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ABSTRACT 


This work used hyperspectral imagery to derive shallow water depth estimates. A 
technique to classify substrates and estimate reflectance values for the substrate types is 
the major contributions of this work. This was accomplished by masking different 
bottom types based on spectra, effects that were not included in previous methods. 
HYDICE data was taken over Lake Tahoe on June 22, 1995. The high altitude of the lake 
provided a low aerosol content within the atmosphere. This allowed for relatively 
straight forward atmospheric corrections. This was substantially easier than in an oceanic 
environment. The atmospheric radiative transfer code MODTRAN3.0 was used to model 
the atmospheric conditions at the time of the experiment. The radiative transfer code 
HYDROLIGHT3.5 was used to model the attenuation coefficients of the relatively clear 
water of the lake. Minimal river input and low chlorophyll concentrations made it 
simpler to determine these values. Making use of the full spectral content of data within 
the optical range, multiple substrates were differentiated and masked off. This allowed 
for an estimation on wet substrate reflectance and a straight forward calculation of bottom 


depth. 
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I. INTRODUCTION 


A basic military need in littoral warfare is an accurate knowledge of near-shore 
bathymetry. This knowledge is necessary for special forces and other combatants prior to 
landing activities, and for marine forces traversing the coastal zone. Such “metoc” 
information is, of course, just one element of the intelligence information needed to plan 
a landing, with other elements including a knowledge of beach trafficability, and shore 
defenses, including mines and obstacles. The work described here addresses how 
bathymetric information can be obtained from (visible) spectral imagery. 

Due to the complex and constantly varying nature of the interaction of 
electromagnetic radiation with water, it’s best to begin the analysis of a new technique in 
a relatively benign environment. Once satisfactory results have been obtained for the 
mode] situation, one can then begin to understand the interaction within the tumultuous 
near-coastal regions of the ocean. Measurement taken with the hyperspectral imager 
HYDICE over Lake Tahoe on June 5", 1995 provided an ideal basis to begin determining 
depth from hyperspectral data. 

As with any measurement of spectral imagery, the data received at the sensor must 
be unmixed with the noise inherent within the medium through which it has traversed. 
For the case of measurements Over water, this noise will include effects due to the 
atmosphere as well as the water column, both of which are extremely dynamic, changing 
with time and geographical position. MODTRAN3.0 is a proven radiative transfer model 
that has been developed over the past two decades and will be shown to provide a 
sufficient model for the Lake Tahoe atmosphere. In addition, the radiative transfer 
model HYDROLIGHT, developed by Curtis Mobley, will be used to determine the 
behavior of the water, or specifically, the wavelength dependent attenuation coefficients. 

This thesis will take previous depth derivation algorithms and build on them to 
take advantage of the wealth of information available through hyperspectral imagery. It 
will conclude by presenting a relatively accurate depth contour of a portion of Lake Tahoe 


called Secret Harbor. [t will begin with a brief presentation of the history of bathymetry 


measurements in Chapter II followed by a discussion of the basic principles needed to 
understand radiative transfer and how light interacts with water in Chapter III. Chapter 
IV will then describe the conditions of the Lake near the time of the measurements and 
how those measurements were taken. Initial observation, analysis and comparison to 
previous algorithms is presented in Chapter V, followed by a complete discussion, in 
Chapter VI, of how to take advantage of the information content within the hyperspectral 
data as it applies to the algorithm. Finally, Chapter VII will present a discussion of the 


results and conclusions drawn from the modeling technique used throughout the thesis. 


Il. BATHYMETRY 


The mapping of the Earth’s oceans dates back to ancient Babylon and times when 
maps were constructed with chisel and rock instead of paper and pencil, or computer, 


Figure 2.1. 


or yo oe a PEP ay acers® 
“te at & nad a * Yeon 
‘ie « Peter, bate be finer 
a0 EO et eg oe FOF S| 
er. yew ge aw. R 
one OE Se ‘er ee 4 "3 , ma. 


: ; ta? * Sey y 
» 4 b a < mee eee oad, idinatt “ale neat 4 
gpePrenaat tah belt ie Sen te 


Ake <7 a ‘ 
fe Sacco Mil es ne emetiie ‘ete 
ete. a 1 a 


6 
tf 

an 
aw 


pr} 
» 





Figure 2.1. Ancient Babylonian 
map depicting Babylon surrounded 
by ocean. Gaskell (1964). 


Figure 2.1 shows an ancient Babylonian map that depicts Babylon surrounded by water, 
somewhat as a castle is surrounded by a moat. This map and those similar to it were 
based on facts they could observe at the time. It wasn’t until Greek mariners and others 
like them bravely and cautiously set out to sea that these ancient ideas on what the oceans 


were like, began to be disproved. 


A. WEIGHTED LINE SOUNDINGS 


One of the first scientific ways in which early mariners could make measurements 


of the ocean depth was with a weighted line, Figure 2.2. 
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Figure 2.2. Depiction of early sounding 
measurements. Gaskell (1964) 


This was an arduous and time consuming method. It often resulted in mediocre depth 
information at best, but until recent times was the only method in use. 

Depth measurements are limited to how much line can be tethered from the 
measuring vessel. In ancient times this meant that measurements were limited to near 
coastal regions. As capabilities of the vessels grew, deeper measurements spanning a 
much larger area of the ocean were possible. As with any measuring instrument, the 
quality of the information produced is a function of the instrument’s resolution. In the 
case of sounding measurements the resolution is, among other things, dependent on the 
quantity of the measurements, how far apart they are made and the ability of the 
measuring vessel to establish an accurate geographical position. In very deep water, as 1s 
normally the case in the open ocean, it sometimes takes several hours to lower and raise 
the sounding dredge. This makes it very difficult to take many closely spaced 
measurements while also maintaining an accurate position. Credit must be given, 


however to the crew of the British ship H.M.S. Challenger. Over the course of 


4 


Challenger’s three-year expedition, the crew made a total of over two hundred soundings 
providing the first look at the relative transoceanic depth. The course taken by the 
H.M.S. Challenger is depicted in Figure 2.3 to provide the reader an idea of the scope of 


the effort put forth by her crew. 
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Figure 2.3. Route taken by the H.M.S. Challenger during it’s three year 
expedition to make transoceanic oceanographic measurements. Gaskell (1964). 


It is interesting to note that in very near coastal water it is more accurate to use a sounding 


pole than a weighted line. 


B. SONAR SOUNDINGS 


With the advent of sonar the same measurements that used to take several hours 
could be made in a matter of seconds. The speed of the measurements allows for a much 
higher frequency of measurement along the ship’s path and therefore a much better 


bottom resolution as shown in Figure 2.4. 
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Figure 2.4. Comparison of soundings taken with weighted line (on the left) and 
soundings taken with sonar (on the right). Gaskell (1964). 


Figure 2.4 is a comparison of weighted line soundings and sonar soundings made of the 
same area of the South Atlantic Ocean floor. As Gaskell (1964) points out, only 13 
soundings were made with the weighted line as compared to the 1300 soundings made 


with sonar, resulting in a much more detailed profile. 


Cc. DEPTH MEASUREMENTS WITH LIDAR 


Just as sonar measures depth using acoustics, a Light Detection and Ranging 
(LIDAR) system use electromagnetic radiation to measure return time. LIDAR however, 
makes use of the different properties of air and water to determine the depth. It operates 
by sending a very short laser pulse downward from an airborne platform. Portions of the 
energy are reflected off the ocean surface and part is reflected off of the sea bed. The 
nature of the interaction between electromagnetic radiation and water will be discussed in 
more detail later in this paper. Given a reasonably distinct bottom return, the depth can 
be calculated by taking the difference between the return times of the surface and bottom 


reflections. 





OPTICAL SURVEY 


ACOUSTIC SURVEY 


Figure 2.5. LIDAR measurements and 
Acoustical measurements, Cassidy (1995) 


As reported in Cassidy (1995), Figure 2.5 displays the results from a test of a French 
system, which shows a comparable accuracy between acoustic and optical results. 
Cassidy argues that a LIDAR has an advantage over acoustical methods in that it is fast, 
allows low cost surveys of difficult to reach or spread out coastal areas. In addition, the 
inherent navigational difficulties associated with coastal sonar surveys are avoided. 
However, it must be kept in mind that as light travels through both air and water, it 
experiences propagation losses that will be discussed in later chapters. This effect in fact 


places limitations on where and how a LIDAR system can be used. 


D. ALTIMETER DEPTH MEASUREMENTS 


Satellite based altimeters are capable of making depth measurements on a much 
wider scale than either sonar or LIDAR as can be seen in Figure 2.6. These depth 
measurements are the result of 4.5 years of U.S. Navy Geosat altimeter measurements 
and 2 years of European Remote Sensing Satellite (ERS-1) altimeter measurements. In 
Figure 2.6 green areas have essentially normal depth, areas with yellow-orange-red hues 


are relatively shallower and areas with blue-violet-magenta are increasingly deeper. 
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Figure 2.6. Depth derivation, on a continental scale, from 
altimeter measurement. From Sandwell et al. (1995). 


As reported by NASA (1986), the sea surface has bulges that result from the variation in 


gravity in different regions of the ocean, Figure 2.7. 
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Figure 2.7. Gravitational effects on ocean 


surface from altimeter measurements. 
From Sandwell et al. (1995). 
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As depicted in Figure 2.7, such features as mid-oceanic ridges have a high concentration 
of mass and therefore will have a greater gravitational pull, causing a “pile up” of water 
above them. This accumulation of water can result in a rise of the sea surface as much as 
5 meters. Contrary, areas were trenches exist will have less of a gravitational pull and 
subsequently cause a depression of the sea surface, sometimes as much as 60 meters. 
These variations in the sea surface can then be measured by an altimeter, using 
electromagnetic radiation, much as a sonar would measure depth via acoustics. 

Altimeter measurements have given scientist an excellent view of the large scale 
depth variation within the Earth’s oceans. However with resolutions on the order of 7 
km, altimeter measurements are not suited for near shore bathymetry where depth 


variations over meter distances are needed. 


Ee. PASSIVE OPTICAL METHODS 


The field of remote sensing can be dated back to as early as 1858 when cameras 
were first placed on balloons and used to take large scale photographs. As outlined in 
Elachi (1987), this was soon followed by kites, then pigeons and eventually airplanes in 
1909. Some of the earliest references that could be found with regard to depth derivation 
from remotely sensed data dated back to World War II, (McCurdy (1940) and Anon 
(1945)). 


1. Satellite Spectral Remote Sensing 


Spectral sensors of the type adequate for littoral or clear water bathymetry are 
relatively few, though the number is set to increase rapidly in the near future. The 
satellite sensors appropriate for this kind of work are the traditional earth resources 
systems, LANDSAT, CZCS (Coastal Color Zone Scanner) and SPOT (Satellite Pour 
Observation de la Terra). Making use of the visible operating range of LANDSAT, 


listed in Table 2.1 (along with other operating characteristics), several papers have 


explored the possibilities for bathymetric depth derivations. 


Table 2.1 Landsat Thematic Mapper Spectral Bands. 
Derived from Collins, 1996 


Band Number Spectral Bands (tm) 


A ae 


a 





1.55 - 1.75 (SWIR) 
10.4 - 12.5 (LWIR) 
2.08 - 2.35 (SWIR) 


In particular Lyzenga (1978) outlines a method of mapping water depth with 





multispectral data. Bierwirth (1993), which will be discussed in more detail later, derives 
an algorithm to get at sea-floor reflectance and water depth by unmixing LANDSAT 
imagery. 

Although no references were found to bathymetric applications for SPOT, it 1s 
very capable to returning data very similar to LANDSAT. Table 2.2 list the different 
operating ranges of SPOT. 


Table 2.2. Operating characteristics for SPOT. Information derived from 
Kramer (1992) 





Mode of Operation Band Spectral Range 


Multispectral 





l 500 nm - 590 nm 
2 610 nm - 680 nm 
3 790 nm - 890 nm 


Panchromatic Black and White 510 nm - 730 nm 
Spatial Resolution = 10m 


10 







The CZCS instrument was Jaunched in 1978 onboard the NIMBUS-7 satellite and was 
the 1° multiple channel optical sensor tuned for observing the ocean environment. CZCS 
data was significant in that it proved that such oceanic constituents as chlorophyll and 
phytoplankton could be determined from remote measurements. However, given a 
resolution on the order of 1 km, CZCS did not prove useful for small scale or shallow 


water measurements. 


Pe Airborne Spectral Remote Sensing 


The Visible / Infrared Imaging Spectrometer (AVIRIS), was one of the first 
airborne spectral imagers. It was developed as a result of the need for greater spectral 
resolution than satellite based instruments could provide and the subsequent high data 
volumes. The success of this sensor prompted a push to develop what is now called the 
hyperspectral sensor and resulted in such systems as the Hyperspectral Digital Imagery 
Collection Experiment sensor (HYDICE) and Advanced Airborne Hyperspectral Imaging 
Spectroradiometer (AAHIS). There are many other instruments currently in operation 
and under development that will not be discussed further here. Hyperspectral systems 
will shortly be included in satellite payloads; the NASA/TRW Lewis satellite is 
anticipated to be the first such, in 1997. 

Initial results from several experiments conducted with hyperspectral sensors have 
been very exciting and have resulted in high quality images as shown in Figure 2.8. 
These data were taken on October 2" over an area of coral reef at Kaneohe Bay, Hawaii. 
The scene was taken by the AAHIJS instrument, operated by SETS Technology, 
Incorporated. AAHIS was the primary instrument flown in the Island Radiance 
experiment conducted by the Hyperspectral MASINT Support to Military Operations 
(HYMSMO) office in October, 1995 staged at Kaneohe Bay, Hawaii. Coincidentally, the 
figure illustrates a number of the problems in the remote sensing area. There is a 
substantial amount of sun glint (small white spots). The substantial color variations 


reflect the variety of bottom types (coral, sand, etc.), as well as water depth. 
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Figure 2.8. Three color image of run 2oct_rll, taken at Island 
radiance. Red - band 50 (705 nm), Green - band 25 (567 nm), Blue 
- band 1 (435 nm). Derived from data provided by HYMSMO. 


Exploitation of the data for water depth was one of the primary goals of the experiment. 
These data offer fair possibilities, but aircraft motion makes geo-registration of the data 
difficult. 

Several experiments have been flown over Lake Tahoe resulting in excellent data. 
Hamilton et al. (1993) applies an empirical model to one of these data sets in an attempt 
to derive depth information. The model used is based on a multiple regression of 
measured parameters, and requires apriori depth information; it will be discussed in 
greater detail later in this thesis. Table 2.3, gives the spectral operating ranges of both the 
HYDICE and the AVIRIS instruments. 
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Table 2.3. Spectral Band Characteristics of AVIRIS and HYDICE. 
Derived from Collins, 1996 


Spectral Range (tum) | Number of Spectral 
Bands 


Kappus et al. (1996) look at Lake Tahoe data taken on June 22°", 1995. They do 










not explore depth derivations, however an initial analysis of the quality and usefulness of 
HYDICE data in determining water radiance parameters is provided. Figure 2.9 from 
Kappus et al. shows that the radiance values determined from HYDICE measurements 


agree closely with the ground truth measurements as well as the modeled values. 


Water Leaving Radiance 





1 
S solid: HYDICE 
c 8 dotted: qround truth 
3 dashed: HYDROLIGHT 
Ps 
o 6 
an 
4 
= 
= 2 
0 = | - a 
0.40 O50 0.50 0.70 


Figure 2.9. Comparison of Remotely sensed HYDICE data to 
that of measured and modeled data. From Kappus et al. (1996) 


As will be shown later, an accurate calculation of the water leaving radiance is one of the 


most important steps in extracting bathymetry. 


3. Recent Developments 


The quality of measurements taken by CZCS prompted the development of 
follow-on instruments such as SeaWIFS to be carried on SeaStar and the Ocean Color 


and Temperature Scanner (OCTS) onboard the Advanced Earth Observing Satellite 
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(ADEOS). ADEOS, considered the follow on to CZCS, was launched in August 1996 
and is dedicated to Earth environmental research. As described by EROC (1996), the 
OCTS sensor will be utilized to observe the ocean environment. Taking advantage of 12 
bands covering the visible and thermal infrared regions, it measures spectral reflectance 
of dissolved substances, phytoplankton and sea _ surface temperature. These 
measurements will be crucial in helping researchers come to a more complete 
understanding of the particulate distribution within water. Understanding this distribution 
better, is a necessary step in deriving shallow water bathymetry. SeaWIFS is expected to 


gather similar information and is expected to be launched in 1997. 
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Ml. OPTICAL MEASUREMENTS 


The taking of optical measurements requires an understanding of (and models for) 
a wide range of optical processes. Atmospheric transmittance and absorption, surface 
reflectance at the ocean surface, and the volumetric scattering all play important roles. In 
addition, when analyzing measurements over shallow waters, reflection off the substrate 
will play an important role as well. In the sections that follow, the optical elements 


needed for this study are presented. 


A. GEOMETRICAL RADIOMETRY 


‘Spectral radiance is the fundamental radiometric quantity of interest in 
hydrologic optics.’, Mobley (1994). It gives a foundation from where all other 
radiometric quantities can be derived, and provides full description of the structure of the 
light field, including the spatial (xX), temporal (t), directional (E ), and wavelength (A) 
dependence. This is in contrast to the irradiance quantities which are measured over all 
directions, and therefore contain no directional dependence. Irradiance describes the 


target illumination while radiance defines instrument measurements. 


L; Radiance 


Equation [3.1] describes the quantities which comprise radiance. AQ is a measure 
of the radiant energy, within the solid angle AQ, that enters a sensor and 1s incident upon 


a detector element of area AA within a time At and over a wavelength band AA . 


_~ = AQ =o | -1 
Hee Ss —— Cw | 34 
i arena 6 mt 5m) ce! 


I5 


pa Irradiance 


In contrast to radiance, when measuring or working with units of irradiance, the 
angular dependence on the amount of radiant energy 1s removed, and the equation is 
reduced to radiant energy per unit time, per unit area, per unit wavelength as in Equation 


[3.2], 


IC = a (Wm? 
Ath\AAX 


nm’). (3.2) 

However, the detectors of interest only receive photons from within a particular 
hemisphere, thus leading to a hemispherical dependence on irradiance measurements. 
While this is a sensor limitation, by rotating the sensor 180°, radiation measurements can 
be made from both hemispheres. For most environmental applications, sensors that 
measure irradiance are positioned straight up to obtain readings of the sky energy - the 
downwelling irradiance, and then straight down to obtain a measure of energy emitted 


and reflected from the Earth’s surface - the upwelling irradiance. 


5. Reflectance 


Two quantities that will be of use are the spectral irradiance reflectance R(z;A) 
and the spectral remote-sensing reflectance R,,(8,0;A), defined as Equations [3.3] and 


[3.4] respectively. 


ENG) 
R(z;i) = ———_., 3.3 
(Z;A) E, (zsh) (3.3) 
R,,(8,0;A) = oe (sr). (3.4) 
| E, (Z=a;h) 


Where E, and Eg in Equation [3.3] are the spectral upwelling and downwelling plane 


irradiance, and R(z;A) is evaluated just below the surface of the water. In Equation [3.4] 
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Ly is referred to as the water leaving radiance and Ey is now evaluated above the surface 
of the water, so that R,, is a measure of the amount of downwelling light that has returned 


through the water surface for detection. 


4, Radiance Invariance 


The radiance invariance law is an important consequence of the measurement. 
Simply stated, “Radiance is distinguished by the property that it does not change along a 
photon path in a vacuum.’, Mobley (1994). This can be illustrated by a geometric 


example, Figure 3.1, showing two different viewpoints of the same system. 


(a) 





Figure 3.1. Radiance Invariance 


In (a) the radiance quotient can be described as ®, / ApQ49, where ®, is the radiant power 
from the surface S,, incident on the collection surface, Ag. The solid angle subtended by 
S, at Ao iS 4249 and distance between the emitting surface and the collector is r. Conversely 
in (b), the radiance 1s described by Dp /A,Q,, where now the radiant power ®p originates 
from a point at the surface S, of variable area A,, and travels within a bundle confined by 


the solid angle 0, to the collector’s surface. In either viewpoint the radiant power 


|e, 


incident on the collectors surface remains unchanged, Mp = ®,. From the definition of 


solid angle, Q = A / i Equation [3.5] follows. 


A,A 
r 


cr 


It then follows from the definition of radiance, that 


L, = ®, [Ap Qo = Do /A,Q, = Lo, (3.6.a) 
thus 


L=L,. (3.6.b) 


In other words, the distance between the source of emission and the collector does not 
change the amount of radiation that arrives at the detector. The relations shown in 
Equation [3.6.a] and [3.6.b] holds as long as the radiation travels within a vacuum. If not 
a vacuum, the medium through which the radiation travels determines how much of the 
emitted signal will be attenuated in the journey to the sensor. With this in mind, models 


can be developed to separate real signals from noise inherent to a particular medium. 


B. LIGHT AND HOW IT INTERACTS WITH WATER 


As light travels through a medium, it will interact in such a way as_ to change the 
characteristics of that light field. Whether these transformations are minor, or extremely 
significant, is dependent on the nature of the medium. In particular the two mediums that 
this paper will be interested in are air and water. The atmosphere, although very dynamic 
and constantly changing, 1s fairly well understood, and several models have been 
developed in the past decades that predict light propagation within it. A brief discussion 
of this interaction and the associated model ‘MODTRAN3.5’ 1s presented in section III.C. 


However, for a more detailed discussion of the subject, the reader is referred to Robinson 
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(1985), or Stewart (1985). However, water 1s a denser medium which contains more 
suspended material in much greater concentration that air. In addition, these 
concentrations change rapidly over very small spatial dimensions making water a very 
difficult medium to model. To understand this interaction, one must first understand how 
the properties of a body of water relate to a light field. Following the reasoning of 
Mobley (1994), the different properties of water can be divided into essentially two 
categories; the first being those properties that depend upon the medium itself, defined as 
inherent optical properties (IOP’s). The second category is composed of those properties 
that depend upon both the medium itself and the directional structure of the light field. 


This second category 1s defined as apparent optical properties (AOP’s). 


1. Inherent Optical Properties 


IOP’s can be better understood by first visualizing how light interacts with a small 


volume of water AV and thickness Ar, Figure 3.2. 
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Figure 3.2. Geometry used to define inherent optical 
properties. From Mobley (1994). 


Using the notation of Mobley (1994), ®,(A) is the incident radiant power of a narrow 
collimated beam of monochromatic light, ®,(A) is the radiant power absorbed by a 


column of water, ®,(A) is a measure of the radiant power that is transmitted through the 
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same column of water, ®,(A) is the radiant power that is scattered by the column of water 
and y is the scattering angle. Summing the different terms in accordance with the 


conservation of energy gives Equation [3.7], 
D(A) = D(A) + BA) + BA). (37) 


From this relation such properties as the spectral absorptance coefficient, a(A), the 
spectral scattering coefficient, b(A), and the spectral beam attenuation coefficient, c(A), 


can be defined. 


a. Spectral Absorptance 


The spectral absorptance is defined as the fraction of incident power 


absorbed within AV, Equation [3.8]. 


_ ®,(A) 


A(r , 
a @ (A) 


(3.8) 


Then by taking the limit of A(A) divided by the length of the water column Ar Equation 
[3.9], 
. AA 


Ar-0 Ar 





with the spectral absorption coefficient a(A) having units of m”. 


b. Spectral Scatterance 


The spectral scatterance is similarly defined as the fraction of the incident 
power that is scattered out of the beam as it passes through the column of water of length 
Ar, Equation [3.10], 


® (A) 


B(A)= 
Y ® (A) 


(3.10) 
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and the spectral scattering coefficient b(A) is defined as Equation [3.11], 





y meee 
b(X) = ]im 9, ; eae: 
Ar-0 r 
C: Spectral Transmittance 


The spectral transmittance, T(A), is given as the ratio of transmitted power 


to incident power as 1n Equation [3.12], 


@, (A) 


T(A) = ; 
ae D(A) 


(Ser2) 


T(A) is a measure of the amount of radiative power that passes through a water column. 


d. Other Significant Quantities 


Several other IOP’s are derived from these 3 quantities. The first 1s simply 
defined as the sum of the spectral absorption and scattering coefficients and is called the 


spectral beam attenuation coefficient Equation [3.13], 
c(A) = a(A) + b(A). Gey 


The beam attenuation coefficient, in turn leads to another important quantity called the 
optical depth, defined as a measure of the attenuation of energy due to both absorption 


and scattering, and given by Equation [3.14], 


Zz 


C= le@ jac? (3.14) 


0 
Where the beam attenuation coefficient c(z) has been expressed as a function of 


geometric depth z. 
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One final quantity of note is called the spectral absorbance — (note not the 
absorptance). This term is more commonly referred to as the optical density and is given 


by Equation [3.15], 


®; (A) 


DO) alone 
1) = 10810 Fa) +O.) 


=—log,,{1- A(A)]. (3.15) 


De Water Constituents 


Knowing the IOP’s is a very important step in being able to model how a light 
field will interact with a body of water. However, these properties depend not only on the 
water itself, but also on the various constituents within the water. It is therefore 
important to be concerned with the various constituents that make up both fresh and sea 
water. The main obvious difference between the two 1s the fact that sea water contains 
various amounts of dissolved salt. Although these salts do not have significant effect on 
absorption in the wavebands of interest, namely the visible portion, they do increase the 
scattering above that of fresh water by approximately 30%. Table 3.1, derived from 
information in Mobley (1994), lists several of the constituents that may be found in both 
types of waters, and gives a brief explanation of each. Particulate matter can, in general, 
be divided into two separate categories based on origin: biological and inorganic sources. 
Those particles that are of biologic origin include bacteria, phytoplankton, zooplankton 
and organic detritus (particulate matter left after the death of an organism and organic 
waste). Inorganic particles enter the water as a result of the erosion of terrestrial rocks or 


soil. 


pays 


Table 3.1. Types of water constituents. 


Colloids Contribute significantly to back scattering 


a Contributes significantly to _ particulate 
backscatter, 

| oe Primarily responsible for determining optical 

roperties of most ocean waters. 









Organic Detritus | Primary backscattering component in_ the 
ocean 
Zooplankton Very small living animals 






3. Summing the Different Inherent Optical Properties 


As described in the last section, water contains many different types of particulate 
matter. Since each of these will interact with a field of light in a different manner, the 
inherent optical properties will change as a function of the distribution of particles within 
a body of water. The water, being a very dynamic entity, also causes the distribution of 
particles to be very dynamic, and therefore difficult to exactly predict. In particular, it’s 
the sum of the effects that is of interest. By knowing the general absorption and 
scattering for different particulate matter, the effects can be summed to develop a feel for 
how the entire body of water will interact with the light field. AOP’s can be generally 
described as a derivative of IOP’s that are dependent on both the nature of the medium 


and the directional structure of the ambient light field. 


4. Absorption in Water 


When discussing the absorption of light in water, most all of the above mentioned 


particulate matter play a role, and need to be modeled. The total absorption coefficient 


ies 


will be the sum of all the different particulate matter coefficients, as well as the inherent 


absorption due to the electromagnetic (EM) properties of pure water. The models 


presented below are taken from Mobley (1994). 


a. Absorption in Pure Water 


For a more complete understanding of EM properties as they relate to 
substances with different index of refraction, the reader is referred to Klein et al (1986). 
For the purpose of this text it 1s assumed that the reader has sufficient background in 
optics to understand the basic principles of plane wave propagation. To begin, a 
relationship between the absorption coefficient a(A) and the complex index of refraction 
k(X) (Also called the Electrodynamic absorption coefficient) is defined by Equation 
[216] 


- Ant k(A) 


a(n) 7 


(3.16) 


Where A is the in vacuo wavelength. Figure 3.3, is a representation of how k(A) varies 


with wavelength. 
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Figure 3.3. Complex (left) and Real Index (right) of 
refraction in Pure Water verses wavelength. From 
Mobley (1994). 
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In addition, Mobley (1994) defines n(A) to be the real part of the complex index of 
refraction m (where m = n - ik). The feature of interest is the nine order of magnitude 
decrease in k(A) as it passes through the near ultraviolet into the visible, and then rises 
sharply again as it approaches the near infrared. This characteristic is directly related to 
the spectral absorption in pure water through Equation [3.16], and is displayed 


graphically in Figure 3.4. 
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Figure 3.4. Absorption coefficient of pure water (solid 
line) and pure sea water (dotted line) plotted verse 
wavelength. From Mobley (1994). 


It 1s this characteristic of pure water and pure sea water that makes shallow water 
bathymetry possible. However it 1s also this characteristic that restricts bathymetric 
measurements to the visible portion of the spectrum. In fact, Figure 3.5 shows that 
shallow water bathymetry is further restricted to the blue - green portion of the visible 


Spectrum. 


25 


0.05 





absorption coefficient a (m“') 
(,W) q Jualoyjaoo Bulsayeos 


0.00 
700 800 


wavelength A (nm) 


Figure 3.5. Absorption coefficient a(A) (solid line) and 
scattering coefficient b(A) (dotted line) of pure water 
plotted verse wavelength. From Mobley (1994). 


Figure 3.5 clearly shows a sharp increase in absorption outside the blue - green portion of 


the spectrum. 


b. Absorption Due to Dissolved Organic Matter 


Dissolved organic matter, which is commonly referred to as yellow matter, 
CDOM or gelbstoff, is well understood. Absorption by yellow matter can be modeled 
fairly closely by Equation [3.17]: 
OR) Sen (Ui) (3.17) 


Thus, by knowing the initial absorption a,(A,) at some characteristic wavelength, 


usually taken to be Ap = 440 nm, the absorption at all other wavelengths can be modeled. 


c Absorption Due to Phytoplankton and Organic Detritus 


Photosynthetic pigments of various types are the major contributors to 
absorption by phytoplankton, of which chlorophyll is known to be the strongest 
contributor. Chlorophyll is common to all photosynthetic plants, and therefore as would 
be expected is a strong absorber of visible light. This effect is particularly strong within 
the green portion of the visible spectrum. Results form the analysis of several different 


types of phytoplankton are plotted in Figure 3.6. 
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Figure 3.6. Total absorption coefficient for selected 
chlorophyll concentrations C. From Mobley (1994). 


It should be noted that each plot of the absorption coefficient takes on a maximum in the 
blue portion of the spectrum at 440 nm and in the red at 675 nm, while it takes on a 
minimum around 600 nm in the green portion of the spectrum. This effect will have 
significant consequences when trying to derive bathymetric information in areas with a 


high chlorophyl!! concentration. 


d. Contributions From Sediment 


Absorption due to inorganic material, although possibly just as significant 
as that of organically derived particles, is not well understood. Much of the research in 
the field of bio-optics has been directed toward understanding the growth of biological 
constituents in different areas of the ocean. As a result, algorithms have been derived that 
help to model the optical interaction of light in waters with varying concentrations of 
organic material. However, comparatively speaking little effort has been put forth in 
trying to better understand the role of sediment in water. However, for the purpose of 


clear fresh water, the effects of sediment load will be assumed to be minimal. 
e. Deriving a Model for Total Absorption 
Several models have been developed that lead to a description of the total 


spectral absorption coefficient for a given water type, each of which in some way or 


pa | 


another will depend on the definitions above. One algorithm, in particularly, will be 
mentioned at this point, due to the consequences it will have on running the radiative 
transfer model Hydrolight3.0, to be discussed later. The algorithm presented in Equation 
[3.18] attempts to model the total absorption a(A), due to all the above terms, and 


express it in one complete formula. 


a(h) = (a, (A) + 0.06a,*"(A)CO% (14+ 022°) — (3.18) 


Sy Scattering in Water 


Scattering can be defined as the redirection of energy, where as absorption is the 
removal of energy. EM radiation can be scattered by virtually any of the above 
mentioned constituents of water. The manner in which it 1s scattered is a function of a 
number of different parameters, including particle shape and size, particle index of 
refraction, wavelength of radiation and viewing geometry. Mie theory characterizes the 
different types of scattering by what is called the scattering size parameter y, Equation 
[3.19]. The Mie parameter is simply a ratio of the circumference of a particle to the 


wavelength A of the incident radiation, 





(Cal) 


where r is the radius of the particle. As would be expected, for different values of x the 
manner in which radiation is scattered will be different. Table 3.2 outlines the different 


ways in which radiation will be scattered as a function y. 
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Table 3.2 Types of scattering based on the scattering size parameter Y. 


As with modeling absorption, it is very difficult to sort out the different individual effects 













within different water types. Therefore, several analytical formulas have been developed 
to model the curves which result from measured data. Several of these models have been 


detailed in Mobley (1994), and will not be discussed in detail here. 


c. RADIATIVE TRANSFER 


Prior to understanding a measured signal, you must come to an understanding of 
the history of the signal, i.e. where it came from and what path it took to arrive at the 
sensor. The theory of radiative transfer explains the rules that govern the transfer of 
radiation from one place to another, whether that transfer is through the atmosphere, 
water or orange paint. The processes that apply to atmospheric modeling are summed 
together in what is called the atmospheric Radiative Transfer Equation, and for our 
purposes can be expressed in elementary form. Using similar notation to that of 


Robinson (1985), the radiative transfer equation is expressed as Equation [3.20], 


Lsensor = Lath a7 lt Lsky + T Lwater. G20) 


Where Leensor iS a measure of the total radiation arriving at the sensor, Loan 1s the 
atmospheric path radiance, L,xy is a measure of the radiance reflected off of the surface, T 


is the atmospheric transmittance and Lyater is defined as the water-leaving radiance. 


a) 


Figure 3.7 depicts the general terms involved and gives a rough idea of the different paths 


the photons take to arrive at the sensor. 





Figure 3.7. Radiative Transfer through the 
air and water. 


The path radiance term is a resultant of all the photons that originated outside of the 
sensor field of view (FOV) and have been, for one reason or an other, scattered by the 
atmosphere into the FOV. The path radiance terms include the dotted lines that lead to 
the arrow labeled Lpath. Radiance that is incident on the water within the sensors FOV 
and is subsequently reflected or scattered back toward the sensor, can be considered to 
come from two general terms -- sun glitter and sky glitter. Sun glitter can be described as 
that radiation from the sun that is reflected at the sea surface, directly into the FOV. 
However, sky glitter is scattered by the atmosphere, prior to being reflected off the 
surface into the FOV. The contributions that make up Ly are represented by the dashed 
lines in Figure 3.7. The final contribution to the overall radiance measured by the sensor 


is defined as the water leaving radiance, Ly, which is shown as the solid lines in Figure 
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3.7. Ly is the sum of those photons that actually enter, interact with and then emerge 
again, from the water, within the sensors FOV. 

As briefly mentioned earlier, MODTRANS3.5 its a radiative transfer model that will 
model the atmospheric effects discussed above. MODTRAN3:5 is the latest generation 
of atmospheric modeling programs developed by Phillips Laboratory. BETA restrictions 
were lifted as of December 1996. At the time this thesis was written, a substantial 
reference for MODTRAN3.5 was not available. However Phillips Laboratory cited Berk 
et al. (1996), which is a paper that had been submitted for publication. A complete 
discussion of the parameters used for the modeling of Lake Tahoe will be provided in 


Chapter IV. 


1. Radiative Transfer at the Water 


Once atmospheric effects are understood and modeled, Ly can be derived. The 
next step is to model the radiative transfer process within the water itself. However, this 
is very difficult as the radiation is effected by scattering at the surface of the water, 
absorption and scattering within the water column and scattering and attenuation of the 
bottom material. All of which are extremely complex and constantly varying. 

Mobley (1994) presents a complete and thorough discussion of the process of 
radiative transfer within the water. HYDROLIGHT 3.0, which was developed by Dr. 
Curtis D. Mobley, can be used to model a variety of different aquatic environments based 
on many of the parameters explained above. The more information a user is able input 
into the code concerning the particular makeup of a body of water, the more accurate the 
results will be. A discussion of the inputs used for Lake Tahoe will be given in Chapter 
V. As will be pointed out later, the parameter of interest, within the water, will be the 
wavelength dependent diffuse attenuation coefficient (Ky). Figure 3.8 1s Figure 8 from 
Mobley (1995) and is provided as an example of HYDROLIGHT’s ability to compute 
Ka. 
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Figure 3.8. Example of modeled Kg. 
Figure 8, from Mobley (1995). 


The values for Kg(A) upwelling (dotted line) and Kg(A) downwelling (solid), plotted in 
Figure 3.8, are calculated for pure water which also contains such particulate matter as 


colored dissolved organic matter and chlorophyll-bearing particles. 


D. BATHYMETRY FROM REMOTELY SENSED RADIATION 


is Unmixing Effects Due to Depth and Substrate Reflectance - The 
Bierwirth Algorithm - An Exploitation of LANDSAT Data 


Water leaving radiance is the amount of radiation upwelling from within the water 
column. This upwelling radiation is the sum of radiation incident off of the substrate and 
the radiation which is due to the bulk reflectance of the water. Let Ly be the water 
leaving radiance, L, the radiance of the wet substrate material, as if there were no water 
cover, Lyw the radiance due to the bulk reflectance of the water column (i.e. deep water 
radiance) and, as before, Ty the transmittance within the medium. Following the method 


Sy 


of Jupp (1988), and Bierwirth et al. (1993), the preceding terms can be combined to form 
Equation [3.21], 


Ly = TwLl; at (le ty) Law - G21) 


This expression now takes into account all of the scattering and absorption of the water 
column by combining them into the one term Ty. If Ty 1s normalized so that it varies 
between O and 1, it then becomes a fraction for the amount of radiation the body of water 
allows to pass. Ty then takes the form of 


T =e 7h, (3.22) 


Ww 


Ky 1s the attenuation coefficient and z is the depth. Equations [3.21] and [3.22] are not 
exact, but come very close to modeling the way in which radiative transfer takes place. 
Jupp, therefore, argues that they sufficiently model the radiative process within the water 


studied. Equation [3.23] is the result of combining Equations [3.21] and [3.22], 


lie Testi Tay = Ioahenuet™. (3.23) 


Bierwirth (1993), follows a similar development to that of Jupp (1988) up to this 
point. Bierwirth then seeks to unmix the effects due to reflectance with those due to 
depth, by deriving a substrate reflectance factor for each band processed. He begins, by 


normalizing Equation [3.23] to reflectance values, as in Equation [3.24] 


R,=R,, +(R, —R,,)e- ea) 


Assuming that the deep-water reflectance is small compared to the substrate reflectance, 


Equation [3.24] can be expressed as Equation [3.25], 


R —R,,, = R 


wt 


1 R ef Kaz. i= l, N. (3.25) 


wi St 
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Where &,,’ is the water leaving reflectance, corrected for deep water, the i-subscripts 
specify a wavelength dependence and N is the number of wavelength measured. Notice, 
that for ideal measurements over deep water, Ryj = Raw; and Equation [3.25] will be zero. 
This indicates that all of the reflectance from the substrate is attenuated. 


Solving for depth z in Equation [3.25], gives Equation [3.26], 


_ In(R,;’)_In(Ry) 
2K, -2K, 


=1,N. (3.26) 

For Equation [3.26], the N values of R,; are unknown as is the depth, giving a total of 
N+1 unknowns that need to be sorted out. A unique solution is unlikely. However, by 
taking a linear combination of all the wavelength dependent equations and assuming the 
geometric mean of the substrate reflectance equals one, a solution for the estimated depth 
falls out. This assumption 1s equivalent to letting the second term on the right hand side 


of Equation [3.26] equal zero. The resultant is Equation [3.27]; 


ENR) 
Z=) —— ~OKN (3.27) 


i=] 
Bierwirth et al. (1993), in applying Equation [3.27] to multispectral LANDSAT data 
taken over Hamelin Pool, Shark Bay, Western Australia has been able to produce an 
estimated depth image, Figure 3.9.a, and 3.9.b. Artificial illumination has been applied in 
Figure 3.9.b to get a better idea of the detailed structure that has been derived. For 
comparison, Figure 3.9.c depicts a true depth image of the pool resulting from 
hydrographic survey. Z, the estimated depth, has an error of Az. A primary contribution 
to this depth error is the assumption of an overall bright bottom, which results in the 


inaccurate report of depth over areas of dark substrate 
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3.9.2.) (3.9.b.) 


Figure 3.9. Derived and measured bathymetry for Hamelin 
Pool, Bierwirth et al. (1993). 


In comparing Figure 3.9.a and 3.9.c, Bierwirth et al. point out that the two are linearly 
correlated “reasonably” well. However, he also makes the point that, in regions of dark 


substrate (i.e. near the bottom and in the in the tidal channels near the top), the depth is 


exaggerated. 


ze Empirical Model - Hamilton Algorithm - An Application of AVIRIS 
Data 


A more empirical approach is possible. Hamilton et al. (1993) estimate depth in 


Lake Tahoe using an empirical model of the form. 


= ay + ay(Rys(A1)) = a2(Ryrs(A2)). (3.28) 


Where R,, is the remotely sensed reflectance at a particular wavelength, as defined by 


Equation [3.4] and ao, a, and az are the linear coefficients. To determine these linear 
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coefficients for Lake Tahoe, Hamilton et al. explain that the surface spectra was taken and 
compared along a transit line of varying bottom depth. Two wavelengths were then 
chosen that displayed a large amount of variance in instrument response, A = 490 nm and 
X. = 560 nm. The application of a multiple regression revealed the linear coefficients to 
be do = 34.96, a; = 23.36 and az = 34.64 with a multiple correlation coefficient of 0.96. 
Figure 3.10.a is a contour plot of AVIRIS data taken over a portion of Lake Tahoe from 


Hamilton et al... For comparison, Figure 3.10.b shows the charted depth of the same 


region. 





(32102) (3.10.b) 


Figure 3.10. A comparison of the depth derived from the 
Hamilton algorithm (3.10.a) to the charted depth (3.10.b) 
for a region of Lake Tahoe. Hamilton et al. (1993). 


As Hamilton et al. concludes, the two scenes are not exact, but do agree in some of the 


major features such as the 60 ft. depth curve. 
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IV. MEASUREMENTS AT LAKE TAHOE 


A. MEASUREMENTS AT LAKE TAHOE 


Due to the complexity of the water environment, as described in the last chapter, 
the best place to start deriving bathymetric estimates with data from passive sensors 1s 
from a relatively clear environment. Lake Tahoe, which 1s located on the California- 
Nevada border high within the Sierra-Nevada mountains provided the ideal conditions to 
begin developing a model for shallow water bathymetry. As explained in Kappus et al. 
(1996), the experiment was conducted on June 22, 1995, on a calm, clear day with very 
little atmospheric aerosols present. Lake Tahoe is 1906 m high at a latitude of 39.14°N, 
and a longitude of 120.19°W, and can be considered a fairly homogeneous body due to 
minimal runoff from rivers and low chlorophyll values (less than .2 mg m”). The lake, 
at the time, was extremely clear. The data was taken at approximately 10:05 am to avoid 
sun glint off of the water, with an aircraft flight path of approximately 100°. This flight 
path was chosen to correspond with the azimuth angle of the sun at the time, again trying 


to minimize sun glint off of the waters surface. 


B. INSTRUMENTS 


Two instruments were flown at Lake Tahoe the day the experiment took place; 
Hyperspectral Digital Imagery Collection Experiment (HYDICE) and the Airborne 
Visible infrared Imaging Spectrometer (AVIRIS). For the purpose of this paper, the focus 
of the discussion will be the use of HYDICE. HYDICE utilizes onboard and laboratory 
based calibration to convert measured raw digital numbers to physical units of radiance. 
HYDICE, which is fitted to fly onboard a Convair-580 aircraft, has a swath that varies 
from 0.2 km to 1.15 km, dependent on the altitude of the aircraft. For the experiment 
conducted over Lake Tahoe, the aircraft was flown at an altitude of 2.35 km (5,000 ft) 


above the lake, which resulted in a swath of 0.385 km and a corresponding resolution of 


oy 


1.2 m. Kappus et al. (1996) points out that this altitude was chosen for a number of 
reasons, including minimizing atmospheric effects, flight path alignment, choice of swath 


width and to allow sufficient collection time. 
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V. INITIAL MODEL APPLICATION 


The HYDICE image of Lake Tahoe was obtained from the HYDICE office on a 
distributed HYDICE 1995 Demonstration Tape, along with several other significant 
scenes taken that year, Table 5.1. Table 5.1 list the tape information on the HYDICE 


1995 Demonstration Tape, with scene of interest in bold type. 


Table 5.1. HYDICE 1995 Demonstration Tape Information 





N/A “support” 
950629 “Yuma, AZ” 2 6514’ 
950622 ‘‘Lake Tahoe, NV” 4 14,544’ 


950622 “Cuprite, NV” 23 25,982’ 
950829 ‘“Aberdeen, MD” 31 OS i 





The scene was written to an 8mm tape in the form of pre-processed radiance data, with 
units of (Watts m” sr’ uum’). Initial analysis and display was completed by utilizing the 
Interactive Data Language (IDL) produced by Research Systems, Inc., Boulder, CO. In 
addition, extensive use of The Environment for Visualizing Images (ENVI), which is also 
a product of Research Systems, Inc., and runs in the IDL environment. The next two 
sections will consist of an initial analysis of the data utilizing the Bierwirth et al. (1993) 
algorithm and the Hamilton et al. (1993) algorithm. This initial analysis is done to get a 
general idea of how the data responds to the different algorithms, thereby giving an 


indication of how to proceed with the final analysis. 
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A. APPLICATION OF THE BIERWIRTH METHOD TO LAKE TAHOE 
DATA 


1. Processed HYDICE Data 


Figure 5.1, is a 320 by 320 pixel display of radiance data taken at Secret Harbor 
on the eastern side of Lake Tahoe, and 1s displayed as a RGB (Red - Green - Blue) image, 
with the Red wavelength set to 650 nm, the Green wavelength set to 550 nm and the Blue 
wavelength set to 450 nm. Selecting the wavelengths in this manner, allows the scene to 
be displayed as a simulated true color image to give the reader a qualitative impression of 


the scene. 








Figure 5.1. Raw data image, displayed with 
EN VI 


Within the scene, shown in Figure 5.1, are regions that have been highlighted by white 
rectangles to display areas in which information for initial data analysis were taken. The 
region labeled 1, was taken over an area of relatively deep water, while regions 2 and 3 


were respectively taken over areas of rocky and shallow substrate. These regions provide 
40 





enough contrast in relation to each other to sufficiently test ability of the algorithm 
presented in Bierwirth et al. (1993), to predict depth. Three pixels were chosen, one from 
each highlighted box in Figure 5.1, for analysis. A line plot of the spectrum for each of 
the three data points, Figure 5.2, indicates a distinct difference in the amount of radiance 
data received at the sensor for each of the different pixels, with the lowest radiance values 
coming from the rocky substrate and the highest values coming from the shallow water 


over a Sandy bottom. 


Lake Tahoe — Hydice Data 
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Figure 5.2. HYDICE Spectra for Three Contrasting 
Pixels. 


Notice how the radiance values for each pixel are centered through the blue - green 


portion of the visible spectrum, as would be expected from previous arguments. 
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Z. Atmospheric Contributions 


The first step in deriving bathymetry is to correct the sensor radiance for 
atmospheric effects. MODTRAN3.5 (Beta version 1.0), radiative transfer code was used 
to model the atmosphere at the time of the experiment. Listed in Table 5.2 are several of 


the parameters used for the different runs made with MODTRAN3.5. 


Table 5.2. MODTRAN3.5 parameters used to model Lake Tahoe 


Run Type Path Radiance Sky Radiance 
Cards Used 1, 1A, 2,3,3A1,3A2,4,5 1, 1A, 2,3, 3A1, 3A2, 4,5 


Parameters 


Visibility 
Ground Altitude 
Initial Height 
Final Height 
Tangent Height 
Viewing Angle 
Scattering 
Julian Day 


Latitude 
Longitude 


Greenwich Time 


Frequency Range 


Frequency Step Size 


120 km 
1.905 km 
4.25 km 
1.905 km 

N/A 
180° 
Mie 
Ie 
39.14° 
120.19° 
17.08z 


10,000 - 25,000 cm’! 


15cm’ 
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120 km 
1.905 km 
N/A 
N/A 
1.905 km 
N/A 
Mie 
173 
39.14° 
120.19° 
17.082 
10,000 - 25,000 cm’! 


15 cm’! 





a. Path Radiance 


Inputting the above parameters, MODTRAN3.5 was first run in radiance 


mode. Assuming that the sensor was looking straight down at the lake, the path radiance 


Lpath Was computed. Figure 5.3 is the resultant spectra for the modcled path radiance. 


Madtran — Path Radiance 


Radiance (W/m? ster s2) 
o> 


O 
0.40 O50 0.60 0.70 0.80 0.90 
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Figure 5.3. Modeled Path Radiance for 
Lake Tahoe, June 22, 1995. 


b. Sky Radiance 


1.00 


Next, the sky Radiance Ly, (described in section II.C) was computed, by 


running Modtran3.5 with the sensor located at the surface, looking up toward the sky. 


The resultant spectra 1s plotted as Figure 5.4. 
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Modtran — Sky Radiance 


Radiance (W/m? ster 2) 
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Figure 5.4. Sky Radiance Computed from 
Modtran3.5, Lake Tahoe, June 22, 1995. 


c: Convolving Modtran3.5 Data to Match HYDICE 


HYDICE has 79 measurement bandwidths of variable width between .4 
and 1.0 um while Modtran3.5 computes the radiance for around 1000 bandwidths. It is 
therefore necessary to convolve the modeled spectra to the spectral coverage of the 
HYDICE sensor. Figure 5.5, is a plot of the sum of path and sky radiance before and 
after it has been converted to HYDICE wavelengths. The convolved spectrum is plotted 


as a solid line. 


ao 


Modtran — Path plus Sky 


Radiance (W/m? ster 2) 
oO 





we 


0 


0.4 0.6 0.8 1.0 
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Figure 5.5. Convolved Path and Sky 
Radiance, Lake Tahoe, June 22, 1995. 


To compute the water leaving radiance Lwyater the wavelength dependent atmospheric 
transmittance Tam 1s also needed. Tatm was computed as part of the path radiance run of 
Modtran3.5. Figure 5.6, is a plot of the modeled transmittance before and after it was 


convolved to the wavelengths of HYDICE. 


Transmittance 


Atmospheric Transmittance 





0.4 0.6 0.8 1.0 
Wovelength (i) 


Figure 5.6. Convolved Transmittance, 
Lake Tahoe, June 22, 1995. 
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d. Water Leaving Radiance 


To compute water leaving radiance for the HYDICE spectrum, it is best to 
look back at Equation [3.20]. However, now the radiative transfer equation must be 


expressed as a function of wavelength as in Equation [5.1], 


Lensor(A) = Lpath(A) + Tatm(A) Lsky(A) + Tatm(A) Lwater(A). (5.1) 


Solving for Lwater gives Equation [5.2], 


E 
water ( alata (r) 


(5.2) 


Substituting the modeled radiance values and the measured radiance of the three selected 


data points into Equation [5.2], results in the water leaving radiance curves of Figure 5.7. 


Hydice — Water Leaving Radiance 





Radiance (W/m? ster 12) 


0.40 0.50 0.60 0.70 0.80 0.90 1.00 
Wavelength () 


Figure 5.7. Water Leaving Radiance, as computed 
for three contrasting data point, shallow-water, 
deep-water and rocky terrain. 
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As with Figure 5.2, the radiance from shallow water (plotted as dark circles with a solid 
line) is the highest, followed by the radiance from deep water (plotted as open circles) and 
the radiance from the rocky region (plotted as a solid line). Recall that at higher 
wavelengths electromagnetic radiation 1s almost completely absorbed (i.e. Ly = 0). In 
fact by observing the radiance values in Figure 5.2 it is seen that this effect is true and 


that corrections for atmospheric effects are reasonably accurate. 


é. Normalizing to Reflectance 


At this point it 1s most convenient to work with reflectance instead of 
radiance. Therefore, following Equation [2.4] a model of the downwelling irradiance 
Eq(A) must be computed. This was again accomplished utilizing Modtran3.5 run in 


irradiance mode. After convolving the data to HYDICE wavelengths the irradiance was 


found as a function of wavelength, Figure 5.8. 


Madtran — Downwelling Radiance 
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Figure 5.8. Down-welling Irradiance, 
modeled for Lake Tahoe, June 22, 


The remote sensing reflectance, R,;, is found by dividing the water leaving radiance, 
Lwater(A), by the down-welling irradiance, Eg(A). Figure 5.9 is a plot of R,; for the three 
data points of interest. 
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Remote Sensing Reflectance 
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Figure 5.9. Remote Sensing Reflectance, Lake 
Tahoe, June 22, 1995. 


The remote sensing reflectance, R,, value will in general vary between O and 1, and. Here 


reflectance varies between O and .018. 


3. Depth Derivation 


Several properties of the water must be known in order to use the method of 
Bierwirth (1993) to derive depth from R,. In particular, to apply Equation [2.27] the 
values for the wavelength dependent attenuation coefficients Kg must be known, or at 
least modeled for the water in question. The radiative transfer model HY DROLIGHT3.0 
was used to carry out the model of the Kg values within the water. Table 5.3 presents 


several parameters that were input into HY DROLIGHT3.0. 
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Table 5.3. HYDROLIGHT3.5 parameters used to model 
Lake Tahoe 


Wavelengths HYDICE Central Wavelengths 
Julian Day 173 

















Latitude a9 14— 










Longitude* -120.19° 


Pressure 





17.65 
5.0 
Humidity 10 % 





Aerosols 





Precipitation 8 % 
Visibility 120.0 km 
Average Chlorophyll 0.2 mg m™ 





concentration 


* West Longitude 1s expressed as negative. 


Several of the above listed parameters are “hard wired”, so to speak, within the file 
‘qarealsky.f’ and must be altered to match the particular environmental conditions of 


interest. 


a. HYDROLIGHT, a Radiative Transfer Model 


Inputting the above parameters, HYDROLIGHT was used to determine 
Kq(A). Figure 5.10, is a plot of the modeled attenuation coefficients that have been 


convolved to the wavelengths of HYDICE. 
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Figure 5.10. Attenuation Coefficients 
Ka(A), as modeled by HYDROLIGHT. 


It should be kept in mind, however, that determination of the attenuation coefficients is 
currently one of the most difficult parts of the problem, and a wealth of research is still 
underway to figure out the best way to model these values. For the moment, this paper 
will assume that HYDROLIGHT has sufficiently modeled Kg, and therefore, use the 
above values to continue the analysis of the Bierwirth method. The sensitivity of the 


results for the Kg values will be examined at the end of the next chapter. 


b. Results of Bierwirth 


Substituting the above modeled values for Kg(A) in Equation [3.27] for 


Kgi, where the subscript ‘i’ indicates the wavelength dependence, the estimated depth can 


be derived for each of the three data points, Table 5.4. 
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Table 5.4. Depths, derived from Equation 
[3.27], for 3 separate data points. 


Bierwirth Depth 









Deep Water 






The values given in Table 5.4 are indicative of assuming a constant bottom reflectance 
R;. The results in the data, as expected contain a large Az factor for each depth. The 
relative depth results for areas over sandy substrate are relatively well behaved. However 
by assuming a uniform bottom type, there 1s no way to differentiate between deep water 
and dark rock. As a result, the calculated depth over the rocks is deeper than that of deep 


water, resulting in unsatisfactory results. 


B. APPLICATION OF THE HAMILTON METHOD TO LAKE TAHOE 
DATA 


In applying the method of Hamilton et al. (1993), the remote sensing reflectance 
Rys, aS explained previously, was computed for the entire scene at similar wavelengths to 
those used in Hamilton et al. Using these computed scene values for R,;, Equation [3.28] 
and the values for ao, a; and a2 given in section III.E.2 a contour plot of the Secret Harbor 
region of Lake Tahoe can be generated, Figure 5.11.a and compared to published 


soundings Figure 5.11.b. 
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Figure 5.11. Comparison of contour plot derived from Hamilton et al. 
(1993) algorithm to published soundings of the same region of Lake 
Tahoe. 


The light colors within this scene represent shallow water as the darker colors represent 
deeper water (white is zero). Therefore, for this scene the Hamilton et al. (1993) 
algorithm has computed the depth to be the deepest near the shore and the shallowest 
further off shore (where white is land). Comparing this to published soundings of the 
same area within Lake Tahoe it is seen that these results are incorrect. To correctly apply 
this algorithm to the Secret Harbor data, apriori depth information would be needed to 
derive the applicable correlation coefficients. However, due to the large amount of 
variation within the substrate of the scene, it would be difficult to establish any solid 
correlation coefficients that apply to the entire scene. Therefore depth results would still 


be incorrect. 


VI. DERIVING DEPTH WITH MODELED BOTTOM TYPES 


As aresult of the discussion in the previous sections, 1t can be surmised that a lack 
of information about the bottom type will result in a depth error when attempting to apply 
the Bierwirth et al. (1993) model to a scene. Therefore, to compute accurate depth 
results, substrate reflectance must be sorted out. This 1s a difficult task when only a few 
wavebands of information are available, such as in Landsat or other. multispectral 
instruments. However, by taking advantage of the wealth of information available in a 
hyperspectral data cube, sorting out the bottom types within a scene becomes much more 


feasible. 


A. MASK CONSTRUCTION 


The HYDICE scene of Lake Tahoe contains a large portion of land, as can be 
viewed in the bottom, right hand, corner of Figure 6.1. The radiant information from 
land, however, is considered noise for the purpose of analyzing the water column within 


the scene and therefore should be removed prior to performing any significant analysis. 
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Figure 6.1. Lake Tahoe, 320 x 320 
pixel data scene. (Displayed at 
wavelength 0.541m) 
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When viewing a scene at wavelengths on the order of .7 um and longer, virtually all of 
the radiation that enters the water is immediately absorbed and never makes it’s way back 
out (i.e. appears black), Figure 6.2. Hence, at these longer wavelengths, the land and the 
water can easily by differentiated, and a mask of the land created, Figure 6.3 (land pixels 


have been set to black). 


. a 
EHS. 
ar 





Figure 6.2. Lake Tahoe, 320 x 320 Figure 6.3. Lake Tahoe - Masked 
pixel data scene. (Displayed at data scene. (Displayed at wavelength 
wavelength | um) 550 nm) 


A plot of the correlation within each waveband, Figure 6.4, for the masked scene, 
now shows what one would expect for a scene that contains only water. Figure 6.4 shows 
that a large amount of information (high correlation coefficient) is contained in the visible 
wavebands around the blue-green portion of the spectrum (.4 - .6 um) and very little 
information (low correlation coefficients) is contained within the red portion of the 


spectrum and beyond (> .6 um). 
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Figure 6.4. Correlation between data 
points within each wave band. 


A principal component analysis (PCA) can now be performed on the masked 
scene to extract significant information. A plot of the eigenvalues of each of the principle 
component bands, Figure 6.5, indicates a high degree of variance within the first few 


bands with very little information in bands 5 and beyond. 
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Figure 6.5. Correlation between data points 
within principal component space. 


aS, 


This four order of magnitude decrease in correlation indicates the relative volume of 
information contained in each transformed waveband. Figure 6.6 shows the relative 


weighting of the observed radiance data contributions to three of the first four PC bands. 
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Figure 6.6. Relative weighting of the 
wavelength dependent data for a few 
of the first PC-bands. (Bands 1,3 & 4) 


The first band, eigenvector 1, is simply a weighted average of the data - roughly the total 


radiance. Bands 3 and 4 include differences which highlight different scene elements. 


1. Constructing Mask for Sandy Bottom Areas 
Figure 6.7 is an image of the first PC band which allows easy distinction of 


shallow water sand and shallow water rock. From this distinction, a mask of the shallow 


water sand is constructed, Figure 6.8. 
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Figure 6.7. PC Band | figure, showing Figure 6.8. Masked shallow water 
contrast between shallow water sand sand in Lake Tahoe scene. 

and rock 

In PC band 4, displayed in Figure 6.9, deep water sand is the prominent feature and is 


contrasted from the surrounding rocks. A mask can be constructed and added to that of 


the shallow water for a total sand mask, Figure 6.10. 
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Figure 6.9. PC Band 4, highlighting Figure 6.10. Total sand mask in Lake 
deep water sand and rock contrasts. Tahoe scene. 
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ie Constructing Masks for Dark Areas 


Thus far, it has been assumed that the dark areas, within the water scene, were 
rocks, however, there are two distinct dark areas as characterized by their different 


radiance spectrum, Figure 6.11. 


Radiance (W/m’st yz) 
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Figure 6.11. Contrasting Spectra from 
dark bottom material. 


The solid line represents the spectra from the brighter of the two material. The difference 
in the spectra suggest that, if the two material are both rocks, then they are undoubtedly 
covered by different material. In fact, this is what will be assumed for the remainder of 
the paper and the two materials will simply be characterized as “bright rock” and “dark 
rock”. 

Taking advantage of the large difference in the radiance around 550 nm, a mask 
can be constructed for the two different types of rocks within the scene, Figure 6.12a and 
6.12.b. This is done by removing all the data within the scene except for that of bright 


and dark rocks, by applying the mask for sand. 
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(6.12.a.) (6. 125.) 


Figure 6.12. Mask for Bright and Dark rocks, (6.12.a and 6.12.b 
respectively) within the HYDICE Lake Tahoe scene. 


The white areas in the scene are the resultant mask for each type of rock. 


oF Composite of the Bottom Types 


Table 6.1 is a summation of the threshold values and bands used to determine the 


different masks. 


_table 6.1. Threshold wees used to define different masks. 


Shallow Sand -13000.0 ~ 3800.0 PC band 1 
Deep Sand 2500.0 PC band 4 


Bright Rock 1000.0 10000.0 Water-Sand Masked | ) 


cube - band36 | 
Water-Sand Masked } 


cube ~ band 36 
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As a result of the different mask formation, the Lake Tahoe scene now becomes a 


composite of three different bottom types; sand - bright rock - dark rock, Figure 6.13. 





masks. (blue - sand, yellow - bright 
rock, red - dark rock). 


The construction of masks for each of these types of bottoms allows for an individual 


analysis of each bottom type and then a reconstruction of the scene. 


B. MODELING DEPTH BY INCLUDING SUBSTRATE REFLECTANCE 


Now that the different areas within the scene can be characterized by the 
reflectance of the underlying substrate, the method of Bierwirth et al. (1993), presented in 
Chapter III, can again be applied to the Lake Tahoe data. However, including effects of 
the substrate reflectance, Equation [3.26] will be utilized instead of Equation [3.27]. 
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ip Estimating Substrate Reflectance 


The process of characterizing the different substrates, was essentially an exercise 
in grouping pixels of similar spectra. Once this is done, depth information can be derived 
by taking advantage of the variability within each of these different substrates 
classifications. To proceed with the application of Equation [3.26] an estimate must be 
made for what Bierwirth et al. (1993) refers to as the radiance of wet substrate material 
for no water cover (L, described in Chapter III.). Ideally this spectra would be determined 
at a depth of zero where the substrate is wet but not covered with water. However, in the 
absence of an exact measurement of the spectra from each of the wet substrates, an 


approximation can be made by taking near shore-values. 


a. Rock Substrate 


The areas within the scene that have been characterized as rock, both the 
bright and dark, have smaller radiance values near-shore than in deeper water. This effect 
is due to the dark material of the rocks reflecting little radiation, similar to a blackbody. 
At shallow depths this dark material will resemble a blackbody more closely than when 
the material is covered by a deeper layer of water (i. e. The bulk water reflectance will 
add to the water leaving radiance over dark areas of the scene). Consequently, the darkest 
pixels needed to be determined in order to obtain a characteristic spectra. A histogram of 
the radiance values from the masked scene can be generated (Figure 6.14 and Figure 
6.15) and, from the values within the minimum bin the pixel with minimum spectra can 
be located and averaged. This was done for both the bright rock and the dark rock, by 
making use of the masks discussed in the previous chapter. For the. bright rock 
approximately 5 pixels associated with a brightness level of around 1050 were selected. 
For dark rock approximately 6 pixels associated with a brightness level of around 750 


were selected. 
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Figure 6.14. Histogram plot for determining Figure 6.15. Histogram plot for determining 
shallow water spectra for bright rock within — shallow water spectra for dark rock within 
the masked Lake Tahoe data cube. the masked Lake Tahoe data cube. 


Notice the difference in radiance values between the bright and dark rock, reiterating the 
need to treat the two as different types of substrate. The resulting spectra are presented 


below. 


b. Sandy Substrate 


In contrast to the dark areas in the scene, the area that has been 
characterized as sand has large near-shore radiance values that decrease with increasing 
depth. Therefore, the best estimate of the wet sand spectra is to determine the average of 
several pixels with the highest spectra. This was again accomplished by plotting a 
histogram of the radiance values as in Figure 6.16. For sand approximately 7 pixels with 


brightness values on the order of 3100 were selected and averaged. 
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Figure 6.16. Histogram plot for determining 
shallow water spectra for sand within the 
masked Lake Tahoe data cube. 


O 1000 


c Wet Substrate Reflectance 


The wet substrate radiance’s are plotted in Figures 6.17.a, 6.17.b and 


6.17.c for wet sand, wet bright rock and wet dark rock respectively. 
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Figure 6.17. Average spectra of wet sand (a), bright rock (b) and dark 


rock (c). 
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The solid line in each Figure represents the average of the various dotted spectra. Notice 
that the bright and dark rock have similar values from about .4 - .5 tm, but differ 
dramatically from about .5 - .6 um. 

In all cases the pixels used to determine the substrate reflectance were located 
near the shoreline. This is the expected result as the water will be the shallowest near the 
shore for each substrate type. Although these values may not be exact for wet substrate 
with zero water cover, they are probably as close as you can get without making ground 


truth measurements of each substrate. 


2: Depth Results 


a. Depth by Using Bottom Reflectance Compared to Depth Without 
Using Bottom Reflectance 


Including the results of substrate reflectance the depth of water in the same 
three pixels used in Chapter V can now be calculated and compared to the results of 


Section V.A.3.b. 


Table 6.2. Comparison of depth derived without substrate reflectance and with substrate 
reflectance. 


Depth without R, Depth with R, 





As would be expected, a relative decrease in depth between the deep water and the rocky 
substrate is now observed. Previous results incorrectly showed depth to be larger over the 
rocky substrate than over deep water. In addition, there is no longer an inherent depth 


error Az as was in the previous calculated depth. These results clearly indicate that to 
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correctly derive depth, the bottom reflectance characteristics must be included tn the 


calculation. 


b. Using Substrate Reflectance to Calculate Depth for Entire Scene 


Equation [2.26] can now be used to compute the depth at each pixel within 
the scene. The results of this computation are displayed as Figure 6.18. Figure 6.19 is 
the published charted depth for Secret Harbor and is provided as a comparison to the 


calculated depth. 
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Figure 6.18. Contour plot of derived bathymetry 
(including effects due to bottom reflectance). 
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Figure 6.19. Published charted 
depth. From NOAA (1987). 


The dark box in Figure 6.19 is roughly the same area as in Figure 6.18, however Figure 
6.19 is shown on a smaller scale to get a better idea of the relative depth of the different 
curves. The red depth curve (or the bottom curve) in Figure 6.19 1s around 4 m. and the 
yellow curve (or the top curve) is around 7 m. Similarly in Figure 6.18 the 4 m and 7 m 
curves are highlighted by a white outline. The depth in Figure 6.19 decreases in 


accordance with the scaling bar to the right of the Figure. 


C. RELIABILITY OF ATTENUATION COEFFICIENTS 


Much of the dependence on accurate depth derivations is dependent on how well 
Kg 1s modeled. Without accurate ground truth measurements for chlorophyll 
concentrations, Ka was the most worrisome parameter used in the calculations. For that 
reason, a brief error analysis was completed on how much a variation in chlorophyll will 
effect the depth calculations. The dynamic chlorophyll range for Lake Tahoe was 


reported in Kappus et al. to be .16 mg m” at the surface varying to .26 mg m”™ at 35 m. 
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Therefore, HY DROLIGHT3.0 was run once for chlorophyll = .16 mg m” and once for 
chlorophyll = .26 mg m”™. The depth was then calculated for the entire scene for each set 
of Kg values resulting in two separate depth scenes. The difference in the depth results 


for these two scenes is displayed in Figure 6.20. 
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Figure 6.20. Variations in depth due 
to variations in Kg values. 


As can be seen from the scaling bar to the right, the error averages around 10% or so and 
goes no higher than 50%. This indicates that the assumption of a relatively uniform 
chlorophyll concentration made earlier will result in only a mild error of depth estimation. 
It is pointed out that the highest errors in the depth calculations result in areas that were 
estimated to be bright rock. This leads the author to believe that these areas are in fact 
rocky substrate that is covered with various amounts of algae. The chlorophyll 
concentrations within the algae then give rise to the higher depth estimation error reported 


in Figure 6.20. 
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VII. SUMMARY AND CONCLUSIONS 


High altitude, clear water and low chlorophyll concentrations made Lake Tahoe 
an ideal spot to begin to develop a method for deriving shallow water bathymetry from 
hyperspectral data. The scene of Secret Harbor on the eastern shore of the lake provided 
a clearly varying substrate type and bottom depth that proved to be a good test case for 
deriving bottom depth. However, given the relatively poor published soundings for the 
scene, it is difficult to fully access how close the depth derivations correspond to actual 
depth. The full data set taken by HYDICE on June 22™ includes the western shoreline 
just south of Dollar Point, which has been charted in much greater detail. It would be of 
interest to apply the methods developed in this thesis to this area of the lake to get a better 
estimation of error. 

Atmospheric conditions were modeled by inputting relevant parameters into the 
radiative transfer model MODTRAN3.5. Similar methods to that of Kappus et al. (1996) 
were used to correct for atmospheric effects with similar results. The clear, virtually 
aerosol free conditions, as described in Kappus et al., allowed for a relatively straight- 
forward modeling of the atmosphere above Lake Tahoe the day of the experiment. As a 
result the derivation of water leaving radiance, Ly, was accomplished with excellent 
results. A marine boundary layer will introduce a larger, and more difficult to account 
for, error. 

The clear water of Lake Tahoe was assumed to be relatively homogeneous for the 
purpose of the calculations in this thesis. Ground truth measurements revealed these 
values to vary between .16 mg m” and .26 mg m” resulting in an average error of about 
10% across the scene for that difference in range. The assumption of a relative 
homogenous chlorophyll concentration should therefore result in an average error of less 
than 10% in depthh HYDROLIGHT3.0 was able to adequately model the values for 
attenuation coefficient Kg given the clear water environment. However, as attempts are 


made to model water with much higher concentrations of particulate matter, it will 
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become increasingly more difficult to model, and closer attention will need to be given to 
thorough ground truth measurements. 

As an initial test of the performance of the depth derivation method stated in 
Bierwirth (1993), the radiance spectra was chosen from three separate pixels were 
analyzed. One of these pixels was located in shallow water over a sandy substrate, one in 
deep water over a sandy substrate and one over a rocky substrate in relatively shallow 
water. Calculations resulted in a depth over the three different areas of interest resulted in 
an erroneous report of depth over the dark rocky substrate with respect to the sandy 
substrate. In addition, each calculation resulted in a large offset error in actual depth. For 
these reasons, it was surmised that substrate effects would need to be included to 
correctly calculate depth. This process has not been included in previous depth derivation 
methods. 

To calculate depth based on substrate type, the scene needed to be divided into 
different regions. A principle component analysis resulted in the classification of regions 
of sandy substrate. Once these regions were masked off, the areas of dark and bright rock 
were differentiated based on spectral differences in selected wave bands. This process 
was easily completed due to the wide selection of spectral characteristics available from 
hyperspectral data. A mask for each of the three regions was created. 

No ground truth data were available on the spectral characteristics of the three 
defined substrates. Therefore, a simple program was developed to select and average the 
most significant near shore values for each bottom type. For sand, the maximum spectral 
values were chosen and for both types of rock the minimum spectral values were chosen. 

Depth was determined for each masked region separately using the respective 
values of bottom reflectance. The three scenes were then added to form one composite 
depth contour. Very good agreement was observed between the derived depth contour 
plot and the published contour plot. 

As a result of the analysis of this thesis, it 1s concluded that it 1s possible to derive 
a detailed map of bottom depth from remotely sensed hyperspectral data. This 1s 


accomplished by the fact that bottom types are distinguishable form one another based on 
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variations within the hyperspectral data. By developing similar techniques to survey 
turbid coastal waters, information can be provided to the war fighter concerning a 
targeted coastal landing zone. This can be accomplished with limited risk to human life 


and military equipment. 
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